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Abstract 

For a class of PT-symmetric operators with small random pertur- 
bations, the eigenvalues obey Weyl asymptotics with probability close 
to 1. Consequently, when the principal symbol is non-real, there are 
many non-real eigenvalues. 
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1 Introduction 

PT-symmetry has been proposed as an alternative for self-adjointness 
in quantum physics Thus for instance, if we consider a Schrodinger 

operator on R", 

P = -h^A + V{x), (1.1) 
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the usual assumption of self-adjointness (implying that the potential 
V is real valued) can be replaced by that of PT-symmetry: 

Vov = V, (1.2) 

where v : R" — )■ R*^ is an isometry with v'^ = 1 ^ v. If we introduce 
the parity operator Uyu{x) = u{v{x)) and the time reversal operator 
Tu = u, then this can be written 

[P, UyT] = 0. (1.3) 

Under mild additional technical assumptions it is easy to see that 
the spectrum of a PT-symmetric operator is invariant under reflexion 
in the real axis. However, in order to build PT-symmetric quantum 
physics it seems important that the spectrum be real, so a natural 
mathematical question is then to determine when so is the case. Re- 
sults on reality and non-reality of the spectrum of PT-symmetric op- 
erators can be found in [121 El [3 [2] . 

The purpose of this note is to show that in a probabilistic sense 
"most" non-self-adjoint PT-symmetric operators that are symmetric 
in the sense of (|2.4p . have their eigenvalues distributed according to 
the Weyl law and hence many of their eigenvalues are non-real. As 
a matter of fact, this will be a rather easy adaptation of general re- 
sults on the Weyl asymptotics for non-self-adjoint operators with small 
random perturbations [HI [TOl [TTl HI [131 [m [5] , where the last three ref- 
erences are the onces that we shall use directly. For technical reasons 
we will state our results for elliptic operators on compact manifolds 
but it would be easy to adapt the results of [13] in order to treat 
Schrodinger operators on R". 

The addition of small random perturbations has the effect of de- 
stroying (uniform) analyticity (if the unperturbated operator has ana- 
lytic coefficients). A very interesting question is to give criteria for PT 
symmetric operators with analytic coefficients to have real spectrum. 

The plan of the paper is the following: In Section [2] we treat the 
semi-classical case and in Section |3| we treat the case of large eigen- 
values. 

2 The semi-classical case 

Let X be a compact smooth manifold of dimension n. Let v : X ^ X 
be a smooth involution; = id, with v ^ id. Fix a smooth positive 
density dx on X which is invariant under v and let us take norms 
with respect to dx. Let P be a a differential operator on X of order 
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m >2 with smooth coefficients so that in local coordinates, 

P= ^c.{x;h){hD,r, a^{-;h)eC°°. (2.1) 

\a\<m 

Here < /i ^ 1 is the semi-classical parameter and we assume that 

aa{x;h) - aa{x;0) = 0{h) (2.2) 

locally uniformly and similarly for all its derivatives. We also assume 
for simplicity that aa{x; h) = aa{x) is independent of h when |a| = m. 
Let 

P{x,€)= Y aa(x;0)C'', pm{x,0= ^ aa(x)C" 

|Q|<m |o|=m 

We assume that Pm{x-,i) / on T*X \ 0, so that P is elliptic in the 
classical sense. We also assume that 

Pm{T*X) / C. (2.3) 

Assume that P is symmetric, 

P = rPT =: P\ (2.4) 

and that 

PU = UP* , where Uu{x) = Uyu{x) := u{v{x)), Tu{x) = u{x). 

(2.5) 

This means that P is PT symmetric: 

[[/r,P]=0. (2.6) 

In addition to the PT-symmetry property (j2.6|) . we have assumed in 
()2.4p that P is symmetric. 

Example 2.1 P = —h'^A + V{x) on T"- where is even and '^V 
is odd, V{—x) = V{x). Then P is symmetric and PT-symmetric with 
v{x) = —X. 

Let R be an auxiliary /i-independent positive elliptic second order 
differential operator on X which commutes with U. We also assume 
that R is real, or equivalently that 

[r,i2] = o. (2.7) 

Then R has an orthonormal basis of real eigenfunctions ej such that 
Uej = {—l)^^^^ej where k{j) = 1 or k{j) = —1. We say that ej is 
even in the first case and odd in second case. Put ej = ej when ej is 
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even and ej = iej when ej is odd. Then {ej} is also an orthonormal 
basis and a hnear combination V = '^ctj^j is PT symmetric iff the 
coefficients Oj are real: U{V) = V. 

In order to formulate our result, we shall follow [T3], where we 
treated a situation without any extra symmetry. 

Let (s C be open, simply connected, not entirely contained in 
:= p{T*X). Let V,{t) := vol{{p e R^*^; \p{p) - < t}). For 
K g]0, 1], z £ Cl, we consider the property that 

V,{t) = Oif"), < t < 1. (2.8) 

Since r i— )■ p{x, r^) is a polynomial of degree m in r with non-vanishing 
leading coefficient, we see that ()2.8p holds with k = l/{2m). 

By B^d{0,r) we denote the open ball in R'^ with center and 
radius r. Let q^j be a random potential of the form, 

qui{x) = ^ ak{uj)ek{x), a{uj) = (afe(t^))o<^jfe<L S B^d{0,R), 

0<Atfe<L 

(2.9) 

where pk > are the square roots of the eigenvalues of h^R so that 
hP'Rej. = p\^k- We choose L = L(h), R = R{h) in the interval 



_ « 2 f (2.10) 

1 < ^ < ^^-M^ M>— -K + (-+ 6)M, 

2 2 

for some e €]0, s — s > ^, so by Weyl's law for the large eigenvalues 
of elliptic self-adjoint operators, the dimension D in ()2.9p is of the 
order of magnitude (L/h)'^. We introduce the small parameter 6 = 
ro/i^i+", 0<To<Vh, where 

' — 71 

TVi := M-hsM-F -. (2.11) 

The randomly perturbed PT symmetric operator is 

Ps = P + 5h^'q^=:P + 6Q^. (2.12) 

Here (cf [13]) the exponent A^i has been chosen so that we have uni- 
formly for h <^ 1 and as above: 

Wh^'^qjL^ < 0{l)h-^/^h''^q^\\H'^ < 0(1), 

where is the natural semi-classical Sobolev space discussed in Sec- 
tion 2 of [H] with a norm equivalent to the standard norm in H'^ for 
each fixed h > 0. 
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The random variables aj{uj) will have a joint probability distribu- 
tion 

P{da) = C{h)e'^'-"'''^L{da), (2.13) 
where for some > 0, 

|V«^| = 0{h~^^), (2.14) 

and L(da) is the Lebesgue measure. {C{h) is the normalizing constant, 
assuring that the probability of i?j^D(0, i?) is equal to 1.) 
We also need the parameter 

eo{h) = (/i'^ + /i" In i)(ln - + (In ^-f) (2.15) 
n tq n 

and assume that tq = to(/i) is not too small, so that eo(/i) is small. 
Recall that Q d C is open, simply connected, not entirely contained 
in Ti{p). The main result of this section is: 

Theorem 2.2 Under the assumptions above, /ei F d have smooth 
boundary, let k g]0, 1] be the parameter in i2.y\) . i2.15]) and 

assume that \2. 8\) holds uniformly for z in a neighborhood of dV. Then 
there exists a constant C > such that for > r > 0, ?> Ce^Qi) 
we have with probability 

>1 _e-^ (2.16) 

^/jn+max(n(A/+l),A''4+J\/) 



that: 



l#KP.)nr)-^^voi(p"Hr))|< 

^ + C(r + ln(^)vol {p-\dV + D{^, r)))) 



(2.17) 



Here #{o'{Ps) D T) denotes the number of eigenvalues of Pg in T, 
counted with their algebraic multiplicity. 



In the introduction of [13] there is a discussion about the choice of 
parameters which applies here also: Very roughly, if tq is equivalent 
to some high power of h, then up to some power of ln(l//i), eq is of 
the order of magnitude /i^. Now choose e = h'^~'^o fpj- gome kq £]0,k[. 
When K > 1/2, then the volume in (f2l71) is 0{r^) with /3 = 2k- 1 > 
and more generally we may assume that it is Oirl^) for some (3 > 0. 
Then we choose r to be a suitable power of h and obtain that the right 
hand side in ()2.17p is 0{h'^~^) for some 7 > 0. With these choices of 
the parameters we also see that the probability in (|2.16p is very close 
to 1. 
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Proof of Theorem 12. 2i We just have to make some small modi- 
fications in the proof of the main result in [TJj (which in turn is a 
modification of the proof in [13]) and only mention the points where 
a difference appears. The proof in the two cited papers (see also the 
lecture notes [15]) uses three ingredients: 

1) The construction of a special perturbation of the form 5q^ with 
Quj as in ()2.9p but with a in the complex ball Bqd (0, R) for which 
we have nice lower bounds on the small singular values of Ps in 
()2.12p . see Proposition 7.3 in [13], Proposition 5.1 in |14j . 

2) A complex variable argument in the a variables using the ex- 
istence of the special perturbation in step 1), which permits to 
conclude that we have nice lower bounds on a relative determi- 
nant for Pg — z, with probability close to 1. 

3) Application of a proposition about the number of zeros of holo- 
morphic functions with exponential growth. (See also [16] for an 
improved version of this proposition, not yet fully exploited.) 

In the present situation we want our special perturbation 6qi^{x) 
to be PT-symmetric, that is we want the coefficients a in ()2.9p to be 
real. All the parts of the proofs in step 1 immediately carry over to the 
case of real a except the following result which is the basic ingredient 
in the iterative process leading to the propositions mentioned above: 

Let ei, Cat be an ON family in L^(X) such that 

JV 

|l5^A,e,-||^^. <0(l)||A||civ 
1 

where the constant 0(1) is independent of the family and especially 
of A^. Then there exists 

0<Hj<L 

with 1 1 a 1 1 CO ^ with the parameters as in (|2.10p . such that 

Mhi < C'(l)/i~^AL"+t+^ 
and such that the matrix 

Mq = ij q{x)ej{x)ek{x))dx)i<j^k<N 
and its singular values 

\\Mq\\ =Sl{Mq) > ... > SN{Mq) 
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satisfy 

llMgll < 0{l)Nh-'^, 

Sk{Mg) > hJ'/0{l), for 1 < /c < N/2. (2.19) 

(See (6.23), (7.20), (7.23) in 

Write q = qi + iq2 where qi = X](3^aj)ej, 92 = s° t\ya.i 

qi and q2 are PT-symmetric. The upper bounds on ||g||_H''' and on 
||Mg|| follow from the bound ||q;|| < R and therefore carry over to qj. 
Since Mg = Mq^ + iMq^ we can apply the Ky Fan inequalities ([S]) 
and get 

W N 

< S2k-l{Mq) < Sk{Mq,) + Sk{Mq,), l<k<—. 



Since the singular values are enumerated in decreasing order, it follows 
that for j equal to 1 or 2, we have 

Sk{Mg^)>^^, l<k<^. (2.20) 

this means that step 1 can be carried out and we get a PT symmetric 
operator Ps as in Proposition 5.1 in [Hj, the only slight difference is 
that rather than taking 9 in ]0, l/4[ we have to confine this parameter 
to the smaller interval ]0, l/8[. 

Step 2 now follows follows from Remark 8.3 in [13], where the 
main point is the reality of the coefficients aj while the assumption 
of reality of the basis elements is not necessary, and was made there 
only because we had in mind a real perturbation. 

Step 3 can be carried out without any modifications. □ 



3 Weyl asymptotics for large eigenval- 
ues 

Let be an elliptic differential operator on X of order m > 2 with 
smooth coefficients and with principal symbol Pm{x,i)- In local coor- 
dinates we get, using standard multi-index notation, 

|a|<m |a|=m 

Recall that the ellipticity of P° means that Pm{x,C) / for ^ / 0. 
We assume that 

Prn{T*X) ^ C. (3.2) 
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As before we assume symmetry, 

(P°)* = rP°r, (3.3) 

and that 

P^U = U{P^)*, (3.4) 

with U = Uy as in Section [21 

Let i? be a reference operator as in and around (|2.7p and define ej 
as there. Write 

Rej = (/x0)2e,-, < ^1! < < < ■•■ (3.5) 

so that = hfi^ where fi^ are given after ()2.9p . Our randomly 
perturbed operator is 

P^ = P'+ql{x), (3.6) 
where lv is the random parameter and 

oo 


Here we assume that a^{uj) are independent real Gaussian random 
variables of variance cj| and mean value 0: 



a° ~ AA(0,c7|), (3.8) 



where 



(/i°)-''e-K')^ <a,<(M°)-^ (3.9) 
M= _„^ ,0</j<-, p>n, (3.10) 

where s, p, e are fixed constants such that 

n n n 
- < s < p , 0<e<s . 

Let H^{X) be the standard Sobolev space of order s. As we saw in 
[5] (where the random variables a° were complex valued), G W^i^X) 
almost surely since s < p—"^. Hence G L°° almost surely, implying 
that has purely discrete spectrum. 

Consider the function F{w) = argpm{w) on S*X. For given 
6*0 G S*^ ~ R/(27rZ), iVo G N := N \ {0}, we introduce the prop- 
erty P(^o, A^o): 

E |V^F(ti;)| ^ on {-u; G S*X; F{w) = 9o}. (3.11) 
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Notice that if P{9o,No) holds, then P{9,No) holds for all 9 in some 
neighborhood of ^o- Also notice that if X is connected and X, p are 
analytic and the analytic function F is non constant, then BNq £ N 
such that ^(6*0, A^o) holds for ah 9o. 

We can now state the main result of this section, which is an 
adaptation of the main result of [S]. 

Theorem 3.1 Assume that m > 2. Let < 9i < 92 < 27r and 
assume that P{9i,Nq) and P{92,Nq) hold for some Nq G N. Let 
g € C°°{[9i,92];]0,oo[) and put 

^k,er,o,x = {^e'^^i <9<92,0<r< Xg{9)}. 

Then for every S g]0, \ — P[ there exists C > such that almost surely: 
3C{oj) < oo such that for all A G [l,c>o[; 

The proof actually allows to have almost surely a simultaneous 
conclusion for a whole family of 9i, 92, g. 

Theorem 3.2 Assume that m > 2. Let Q he a compact subset of 
[0, 27r]. Let A'o G N and assume that P{9,Nq) holds uniformly for 
9 £ Q. Let g be a subset of {{g,9i,92); 9j £ 9,6*1 < O2, g £ 
C°°([^i) ^2]; ]0, oo[)} with the property that g and 1/g are uniformly 
bounded in C°°([0i, ^2]; ]0, oo[) when (3,^1,^2) varies in Q. Then for 
every S g]0, ^— /?[ there exists C > such that almost surely: 3C{uj) < 
00 such that for all A G [l,oo[ and all (fif, ^1,^2) S Q, we have the es- 
timate (MM- 

The condition (|3.9|) allows us to choose aj decaying faster than any 
negative power of /i^. Then from the discussion below, it will follow 
that qu){x) is almost surely a smooth function. A rough and somewhat 
intuitive interpretation of Theorem 13.21 is then that for almost every 
PT symmetric elliptic operator of order > 2 with smooth coefficients 
on a compact manifold which satisfies the conditions (|3.2p , ()3.3p , (|3.4p , 
the large eigenvalues distribute according to Weyl's law in sectors with 
limiting directions that satisfy a weak non-degeneracy condition. 

Proof of Theorem 13.11 As already mentioned, the theorem is a 
variant of Theorem 1.1 in [5]. The difference is just that we now use 
real random variables in the perturbation in order to assure the 
PT-symmetry while in [5] they were complex. The proof in [5] used a 
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reduction to the semi-classical case where the main result of jl4j could 
be applied. The proof of Theorem 13. II is an immediate modification of 
that proof, where we replace the main result in [T3] by Theorem 12. 2[ 
The only point where the use of real Gaussian random variables in 
stead of complex ones causes a slight change is the use of (4.10) in [5] 
that was established in [5], where we have to replace the denominator 
2 by 4 in the case of real random variables. That was also proved by 
Bordeaux Montrieux in [3], Proposition 2.5.4. □ 
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